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Abstract 

Clonal interference, competition between multiple co-occurring beneficial mutations, has a major 
role in adaptation of asexual populations. We provide a simple individual based stochastic model of 
clonal interference taking into account a wide variety of competitive interactions which can be found in 
nature. In particular, we relax the classical assumption of transitivity between the different mutations. 

It allows us to predict genetic patterns, such as coexistence of several mutants or emergence of Rock- 
Paper-Scissors cycles, which were not explained by existing models. In addition, we call into questions 
some classical preconceived ideas about fixation time and fixation probability of competing mutations. 
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Introduction 

From the works of the ’great trinity’ of Fisher OHa], Wright [45] and Haldane [22] the questions of 
fixation probability and fixation time of new beneficial mutations have been widely studied. These are 
indeed fundamental questions if we aim at understanding how and how fast a population can adapt to a 
changing environment, the dynamics of genetic diversity, or the long term behaviour of ecological systems. 

The first models of adaptation in asexual populations postulated that the beneficial mutations were 
rare enough for populations to evolve sequentially by rapid fixations of a positively selected mutation 
alternating with long periods with no mutation (see [38] for a review of these models). However, various 
empirical evidence [nisniisi] show that in large asexual populations, several mutations can co-occur, 
which especially can lead to a competition between beneficial mutations. This phenomenon is known as 
clonal interference |21j . and has a major importance for adaptation of asexual populations, such as bacteria 
and other prokaryotes, yeasts and other fungi, or cancers. Consequently in recent years there has been a 
growing interest in developing experimental studies and theoretical models to analyse clonal interference 

[i[iii[ii[n[[iii[Mii3n]. 

Most models investigating how clonal interference might affect the probability and time of fixation of 
beneficial mutations, and thus adaptation, made two important but limiting assumptions: first population 
sizes are constant and independent of the fitness of the individuals, and second fitnesses only depend on 
the type of the mutations and not on the state of the population, i.e. fitnesses are assumed transitive. 
However, as emphasized by Nowak and Sigmund m, there is a reciprocal feedback between adaptation 
and environmental changes because of ecological interactions between individuals, especially competition. 
Mutation with the highest fitness invades the population, but its fitness might depend on the density of 
the other mutations present in the population and might thus change during the course of adaptation. 
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Thus, fitnesses are not necessarily transitive, i.e. selection can be frequency dependent, and phenomena 
such as cyclical dynamics or stable coexistence can occur. Interestingly, non-transitivity of fitnesses has 
been documented empirically in asexual populations [39l [28] . 

In this paper we aim at providing and studying a simple individual based stochastic model taking into 
account the wide variety of competitive interactions which can be found in nature. We show that type 
dependent competitive interactions are able to generate ecological patterns which are observed but not 
explained by conventional models. In particular we relax the classical assumption of transitivity between 


the different mutations (see the discussion in Section 2.1) and call into questions some classical preconceived 
ideas about clonal interference (see Propositions!^ and. To model precisely the interactions between 
individuals we extend the model introduced in [^where the author only considered the occurrence of one 
mutation. The population dynamics, described in Sectionj^ is a multitype birth and death Markov process 
with density-dependent competition. We reflect the carrying capacity of the underlying environment by a 
scaling parameter iti G N and state results in the limit for large K. 

Such an eco-evolutionary approach has been introduced by Metz and coauthors [,821 and has been made 
rigorous in the seminal paper of Fournier and Meleard [5D]. Then it has been developed by Champagnat, 
Meleard and coauthors (see [ziiniis] and references therein) for the haploid asexual case, by Smadi [32] for 
the haploid sexual case, and by Collet, Meleard and Metz m and Coron and coauthors HKll] for the 
diploid sexual case. 

This work is the first one to study the impact of type dependent competition on clonal interference. 
By using couplings with birth and death processes and comparison with deterministic system, we are able 
to provide a complete description of the possible population dynamics when two mutations are in compe¬ 
tition. We show that a large variety of behaviours can be observed, but that in all cases, the dynamics 
consists in an alternation of long stochastic phases (with a duration of order log IF) and short phases 
(with a duration of order 1) which can be approximated by deterministic processes. In the case of Rock- 
Paper-Scissors dynamics, we can even have a non bounded number of such alternations (see Proposition 
[^. The deterministic approximations will be either two-dimensional or three-dimensional Lotka-Volterra 
competitive systems. Unlike the two-dimensional case where the final state is completely determined by 


the signs of the invasion fitnesses (defined in (1.6) and (1.16)), the long time behaviour of a three dimen¬ 
sional competitive Lotka-Volterra system can depend on the initial state. Moreover, the flows of three 
dimensional Lotka-Volterra systems do not necessarily converge to a stable equilibrium but can exhibit 
cyclical behaviours. We will see that the time of appearance of the second mutation is crucial to determine 
the final population state. For example, the dynamics described in Proposition can only be obtained if 
the second mutation occurs when the fraction of the wild type-individuals in the population is small. 

In Sections [T] andwe present the model and the main results of the paper. Sections]^ toare devoted 
to the study of the dynamics of stochastic and deterministic phases. Section l^is dedicated to the proofs 
of the main results. Finally in the Appendices we present technical results |A] ) and provide a complete 
description of the possible population dynamics (]b|. 


1 Model 

We consider an asexual haploid population and focus on one locus. The individuals may carry three 
distinct alleles 0, 1 or 2, and we denote by £ := {0,1, 2} the type space. The population process 

= {N^{t),t>0) = ((Vo^(t),7Vf(t),Vf(t)),t > 0), 

where {t) denotes the number of f-individuals at time t when the carrying capacity is K, is a multitype 
birth and death process. The birth rate of i-individuals is 


h{n) = Piiii, (1.1) 

where /3i > 0 is the individual birth rate of f-individuals and n = (no,ni,n 2 ) G denotes the current 
state of the population. An i-individual can die either from a natural death (rate <5^ > 0), or from type- 
dependent competition: the parameter Cij > 0 models the impact an individual of type j (resp. i) has on 
an individual of type i, where (j,j) G The strength of the competition also depends on the carrying 
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capacity X S N, which is a measure of the total quantity of food or space available. This results in the 
total death rate of individuals carrying the allele i G £: 

(n) = I Oi + —^'^0 + + “^*^2 I rii. (1.2) 

Let us now introduce the notions of invasion, final state and time of invasion that we will use in the 
sequel. We say that a type i G £ invades if there exist h G P,V > 0 and iG G N such that the following 
event holds: 

I{i, K,fi, (3,V) := {y f3\ogK < t < , N^{t)>nK}. (1.3) 

When K is large, it means that i-individuals represent a non negligible fraction of the population during a 
very long time. It is necessary to bound this duration (by as a birth and death process with density 

dependent competition almost surely gets extinct in finite time. To characterize the long time behaviour 
of the 0-, 1-, and 2- population sizes, we introduce the following notion of final state set for n G 
e,l3,V >0 and K gN: 

{nGFSie,K,l3,V)}-.= \ sup \\N^{t)/K - n\\ < e], (1.4) 

'' p log K<t<e''^ '' 

where ||.|| denotes the L^-Norm on Finally, we introduce the time needed for the rescaled population 
process /K to hit a vicinity of h G and stay in the latter during an exponential time. For e,V > 0 
and iF G N: 

Tp^y){n) := inf{s > 0,Vs < t \\N^{t)/K — n\\ < e}. (1.5) 

As a quantity summarizing the advantage or disadvantage a mutant with allele type i has in a j- 
population at equilibrium, we introduce the so-called invasion fitness Sij through 

Sij ■■= (3i- Si - Cijfij, (1.6) 


where the equilibrium density is defined by 


rii 


A - Sj 
Cr,. ■ 


The role of the invasion fitness Sij and the definition of the equilibrium density rii follow from the properties 
of the two-dimensional competitive Lotka-Volterra system: 


= (A - Si - = Zi 


n, 




n-j ' = (A - A - Cj,^n'l^^ - Cj,jn'f’)n'f’, 


Tz)/ 


A A) = 2:^ 


(1.7) 


for z = {zi,Zj) G K+. If we assume 


rii > 0, rij > 0, and Sji < 0 < A 


( 1 . 8 ) 


then rii is the equilibrium size of a monomorphic ^population and the system (1.7) has a unique stable 


equilibrium {ni,0) and two unstable steady states (0,nj) and (0,0). Thanks to Theorem 2.1 p. 456 in 
[I^ we can prove that if Nf'(O) and AfA(O) are of order K, K is large and N^{0) — 0 iov k = £ \ {i,j}, 


the process /K,N^/K) is close to the solution of (1.7) during any finite time interval (see (A.4) 


for a precise statement). The invasion fitness Sij corresponds to the per capita initial growth rate of the 
mutant i when it appears in a monomorphic population of individuals j at their equilibrium size [njK\. 
Hence the dynamics of the type i is very dependent on the properties of the system (1.7). It is proven in 


[7] that under Condition (1.8) one mutant i appearing in a monomorphic _)-population at its equilibrium 
size [njK\ has a positive probability to fix. More precisely, for e > 0, there exist two finite constants /3(e) 
and V (e) such that 


(^nid G FS{s,K,/3{s),V{e)) N^{0) = ei+ [njK\ej'^ = ^ + Ok{£), 


(1.9) 
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where {ei,i G £) is the canonical basis of and Ok{£) is a function of K and e satisfying 

limsup |Oic(£)| < ce, (1-10) 

K—^oo 


for a finite c. Similarly, 


lim P( njCj G FS{e, K, j3{e),V{e)) N (0) = + \ njK\ej 1=1-^ + Ok{£)- 

K—^OO ' ! f-i 


s,.. 


( 1 . 11 ) 


Moreover, the invasion time of the mutant i, as defined in (1.5), satisfies for a finite c 
logiG 


((1 -ce)^< < (1 + ce)^\N^{0) = e. + ln,K\e,) =^+ OK{e). (1.12) 


If we assume on the contrary 

fii > 0, Uj > 0, Sij > 0, Sji > 0 and Ci^iCjj ^ CijCj^i, 


(1.13) 


then the system (1.7) has three unstable steady states {ni,0), {O.rij) and (0,0) and a unique stable 
equilibrium (n|*^, , where 


ij ; "zj 

-(i) _ ~ ^i) ~ — Sj) 


= 


Ci^zCjj - Cz,jCj,z 


and ^i) 

and n C — C C 


(1.14) 


It is proven in that under Condition (1.13) one mutant i has a positive probability to invade then 
coexist with a j-population during a long time. More precisely, for e > 0, there exist two finite constants 
13 (e) and V{£) such that 

P(h«e, +hg\, G FS{£,K,p{e),V{e))\N^{{)) = e, + = ^ + Ox(£), 

P[fije^- G FS{£,K,P{e),V{e)) iV^(O) = + [njK\eA = 1 - ^ +Ok{£)- 

^ 'Pi 

Moreover, the invasion time of the mutant i, as defined in (1.5), satisfies for a finite c 
P(( 


(l-c£)i^ 






\ogK 




iV^(O) = e. + \n,K\e,) = ^ + Ok^s). 


Analogously, the fate of a mutation of type k occurring when the types i and j coexist in a population 
of large carrying capacity K, or the result of the competition between the three types of populations when 
they have all a size of order K is very dependent on the properties of the three-dimensional competitive 
Lotka-Volterra system: 


= (/3i - Si - 




n- '(0) = Zi, 

= iPj - Sj - ^( 0 ) = Zj, 

= iPk -5k- Ck^in'P - Ck^jU^p - Ck,kn^k'')'^^P ^ 


(1.15) 


where z = (zi^Zj^Zk) G Similarly as in (1.6l we can define the invasion fitness of the type in a 
two-type z/j-population, 

(1.16) 


Skij — Skji ■— l^k ^k Ck^iTT'^} Ckjfl) 


where we recall Definition (1.14). It corresponds to the initial per capita growth rate of an individual of 
type k appearing in a large population in which the individuals of type i and j are at their coexisting 
equilibrium, ([n-*'’Arj, [nlj^ATj). 

The class of the three-dimensional competitive Lotka-Volterra sytems has been studied in detail by 
Zeeman and coauthors [4811471146] . They exhibit much more variety than the two-dimensional ones, and 
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we will recall the long time behaviour of their flows in Section 4.2 The complexity of the three-dimensional 
competitive Lotka-Volterra systems entails a broad variety of dynamics in the case of two mutations oc¬ 
curring successively in a population. In particular, depending on the population state when the second 
mutation occurs, the new type can be advantageous or not when rare: for example, we can have S 20 > 0 
and 5'2 io < 0, which implies that if the mutant 2 occurs in a 0-population at equilibrium the 2-population 
size has a positive probability to hit a non-negligible fraction of the population, whereas if it appears in a 
two-type 0/1-population at its coexisting equilibrium it gets extinct immediately. 


In the sequel, we assume that a mutant of type 1 occurs in a population of type 0 (also called wild 
type) at equilibrium at time 0, and that the mutants 1 are beneficial when rare. In other words we assume 


Assumption 1. 


( 0 ) = lnoK\eo + ei and S'lo > 0. 



Figure 1: The three phases of a mutant invasion. The y-axis corresponds to the two type population sizes 
(0 in black, 1 in red), and the x-axis to the time. In this simulation, K = 1000, (/3o,/3i) = (2,3),5i = 
0.5,C,,,=l,(i,j)G{0,l}2. 


Let us describe in few words what happens when there is only the first mutation (we refer to [7] for 
the proofs). We can distinguish three phases in the mutant invasion (see Figure [^: an initial phase in 
which the fraction of mutant-individuals does not exceed a fixed value e > 0 and where the dynamics of 
the wild-type population is nearly undisturbed by the invading type. A second phase where both types 
account for a non-negligible percentage of the population and where the dynamics of the population can 


be well approximated by the system (1.7). And finally a third phase where either the roles of the types 


are interchanged and the wild-type population is near extinction (when (1.8) holds), or the two types stay 
close to their coexisting equilibrium (when (1.13) holds). The duration of the invasion process is of order 
log A with a deterministic phase which only lasts an amount of time of order 1. More precisely, during 
the first phase, when the population 1 has a size smaller than \eK \, this latter can be approximated by 
a supercritical birth and death process, with respective individual birth and death rates 


/3i and 


(5i -I — — [3i — ^lo < Pi- 

K 


Hence from well known properties of supercritical birth and death processes (see [T] for example) we 
know that with a probability close to Sio/Pi the 1-population size hits the value \sK\ in a time of 
order log A/ 5'10 for large A. The second phase can be approximated by the dynamical system (1.7) with 
(bj) = (0,1) and the system takes a time of order 1 to get close to = ([h^Aj, AJ) (if 

5oi > 0) or to = ([e^Aj, [hiAj) (if 5oi < 0). Finally, if 5oi > 0 the populations 0 and 1 stay 

close to their coexisting equilibrium ([hgi^Aj, \ n^QiK\) during a time of order for a positive V and 
if 5oi < 0 the 0-population size is comparable to a subcritical birth and death process, with individual 
birth and death rates 


Pq and 60 + 


—^mK — Pq + |5oi| > Po- 
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Hence it gets extinct almost surely, in a time close to log-ftr/|S'oi|. 


To summarize, for large K, the first and third phases have a duration of order at least log K and the 
second phase has a duration of order 1. Hence if a second mutation appears during the sweep, this occurs 
with high probability during the first or the third phase. We denote by a log K the time of occurrence of 
the second mutation and we distinguish three cases: 

• Either the mutation 2 occurs during the first phase of the mutant 1 invasion and the 1-population 
size hits [sK\ before the 2-population, which corresponds to the Assumptions and: 

Assumption 2. 

5'2o > 0 and 0v(-^;-—) < a < 

VOio *10 


• Or it occurs during the first phase but the invasion fitness S 20 > 0 is large enough for the 2-population 
to hit size [eK\ before the 1-population, which corresponds to the Assumptions and: 


Assumption 3. 


S'20 > 0 and 0 <a < ls2o>Sio (- 

VDlO *20^ 


Or it occurs during the third phase, which corresponds to the Assumptions and: 

Assumption 4. 

1 


a > —— and 
<-’10 


S'oi > 0 
or 

Sqi < 0 and a < 1/S'io + l/|5'oi|- 


In Appendix]^ we give a complete description of the possible population dynamics. We now highlight 
the biologically relevant outcomes. 


2 Results 

2.1 Transitive versus non-transitive fitnesses 

Let us say a word about classical models of clonal interference (see [3 [211 m 1 [H] for instance). The 
population size is constant (or infinite). The fitness Si of an z-individual corresponds to its exponential 
growth rate (Malthusian parameter) and only depends on its type i. Suppose that in a population of type 
0 a beneficial mutation 1 is followed by a second beneficial mutation 2 before the fixation of the type 1 
individuals. In population genetics, ’’beneficial” means that si > 0 and S 2 > 0, which implies that the 
mutant populations have a positive probability to escape genetic drift and constitute a positive fraction 
of the population. Then there are two possibilities: 

1. Either si > S2: then the I-population outcompetes the 0- and 2-populations and the mutation I 
becomes fixed. 

2. Or S 2 > Si: then the 2-population outcompetes the 0- and 1-populations and the mutation 1 becomes 
fixed. 

If a third mutation (individuals of type 3) occurs with fitness satisfying S3 > S2, then not only the type 

3 individuals outcompete the type 2 individuals, but by transitivity of the total order > on K, they also 
outcompete the type 1 and type 0 individuals, and so on. In other words, in population genetics of 
haploid asexuals, the fitnesses are transitive in the sense that if 1 outcompetes 0 and 2 outcompetes 1, 
then necessarily 2 outcompetes 0. 

Such a model is natural when competitive interactions between individuals are simple: in an exper¬ 
iment with only one limiting resource, benehcial mutations often correspond to an increase of resource 
consumption efficiency. But let us imagine an environment with two resources, A and B. A mutant which 
prefers resource A (resp. B) will be favoured in a population of individuals which consume preferentially 
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resource B (resp. A). In this example there is no transitivity. Moreover, experiments on cancer and viral 
cells have shown the long term coexistence of several mutant strains |25| , which cannot be explained by a 
model with transitive fitnesses. 

More generally, it is known m that ecological interactions cause non-transitive phenotypic interactions. 
Hence to model complex interactions, we need an other definition of ’’fitnesses”, to make appear the 
dependency on the population state. This is achieved by the notion of invasion fitnesses (defined in (1.6) 
and (1.16)), which naturally follows from the individual based model that we have presented. 

The main novelty of our approach is to consider type dependent competitive interactions. Indeed, if 
all the competitive interactions {Cij, (i,j) G have a same value C, then for three individual’s types i, 
j and k, the invasion fitnesses satisfy: 


> 0 


Sji = /3j - = /3j - 6j - CiPi - Si)/C = -Sij < 0, 


and 


Sij > 0 and Sjk > 0 


Sik — Sij -f S- 


'jk 


> 0 . 


In other words, it boils down to the case of transitive fitnesses. We can have a more precise result on the 
form of competitions allowing non transitive relations between several mutants. To state this latter we 
introduce the following order for {i,j) G £^: 


i~<j 4 
i=j 


Sij < 0 < Sj i , 

» Sij.Sji > 0 , 


and the notation 


^ • 


( 2 . 1 ) 


Then we can state the following Lemma, which will be proven in Appendix [A| 
Lemma 2.1. Let Ci < C 2 be two positive real numbers and i-,j,k G £. 

1- If 

1 " 


(CiV—) >C 2 and (^AC 2 ) < Ci 


( 2 . 2 ) 


then if for every m n G {*, j, k}, Ci < Cm,n < C 2 , 

i < j and j <k 

2. If 


i ^ k. 




> C2 




<Ci 


(2.3) 


then if for every m n G {*, j, k}, Ci < Cm,n < C 2 , 

i < j and j <k 

3. If 


i < k. 


(--A)- 


< C 2 and ( 




>Ci 


(2.4) 


then there exist some ecological parameters {prmCm,m'rn n G {*,i, fc}) such that for every m 
n G {i,j, fc}, Cl < Cra,n < C 2 and 

* ^ j ^ k o-nd k ^ i. 

Remark 1. By looking at all the possible subcases we can show the following equivalencies: 

CI<C 1 <C 2 <I or 1 < Cl < C 2 < Cl, 

Cl < C| < C 2 < 1 or 1 < Cl < Ci^ < C 2 , 


( 2 . 2 ) 


(2.3) 


(2.4) 


Cl < 1 < C 2 . 


This shows that transitive relations are more likely when the rescaled competitions {Cij,{i,j) G £^) are 
close to each other and both smaller of greater than 1. 
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By allowing dependency of competitive interactions on the individual’s types, we are able to model 
new ecological patterns found in nature and describe them mathematically. In particular, we will show 
that it allows us to model coexistence of several mutant populations, and complex dynamics as Rock- 
Paper-Scissors cycles. 

To present our main results in a simple way, we introduce the following notation, 

p(fe)( ) =p(.|the mutation(s) k is(are) present), k S {{!}, {2}, {1, 2}}, 
where we recall that the mutation 1 occurs at time 0 and the mutation 2 at time a log K. 


2.2 Does clonal interference speed up or slow down invasion? 


In population genetic models of clonal interference, the authors concluded that the presence of other 
mutants slowed down the invasion of the fittest mutant as it had to outcompete some individuals fitter 
than the wild type individuals [33 [Ml HD [HI- In our individual based model, the result of the competition 
between two mutant subpopulations depends on the state of the total population, and the interplay between 
different mutations can take various forms. For some values of the parameters, clonal interference does 


indeed slow down the invasion of a benehcial mutant. Recall Equations (1.91 to (1.12) about the invasion 
of one mutant. Then 


Proposition 1. Let Assumptions^^and^he satisfied and suppose: 

'S'oij*S'o 2 < 0, Si 2 < 0 < S 21 , and S 21 < S 20 , 


(2.5) 


or 

S'oi > 0,iS'o2 < 0, S 12 < 0 < 5'2oi and S '201 < S 2 o- (2-6) 

Then the presence of the mutation Idoes not modify the invasion probability of the mutation 2 and the 
final state set of the population: for e > 0, there exist two finite constants l3{e) and V{e) such that 


p('=)(h 2 e 2 e FSie,K,/3ie),Vis))) = ^+OfiK), for k = {1,2} or {2}. 

P2 

However the presence of the first mutant slows down the invasion of the second one: there exists a positive 
constant c such that for every e > 0; 


- ce)i^^ <Tpg^^^^^^{n 2 e 2 ) - alogK < (l-Fce)^^^) = AOr^s), 


J20 


520 


where 82 ^ > S 20 is defined by 


520 S^\Sio 5 J 52o’ 


where S = S 21 in case (2.5) and S 201 in case (2.6), and 

log if 


P^«((l - < (1 + = 1^(1 - {f) + OKi.y 


) = —I 

fl-—) 

^ P 2 



But due to the non-transitive phenotypic interactions that our model is able to take into account, the 
presence of the first mutant can also speed up the invasion of the second one. More precisely we have the 
following result, which is illustrated in Figure 

Proposition 2. Let Assumptions^^and^be satisfied and suppose: 

'S'oii5'o2 < 0, 5i2 < 0 < 521, and S 21 > S 20 , (2-7) 

and 5201 > 52o- (2-8) 


or 


5oi > 0, 5 o2 < 0, 5i2 < 0 < 5201 










Then the presence of the mutation Idoes not modify the invasion probability of the mutation 2 and the 
final state of the population: for £ > 0, there exist two finite constants j3{e) and V{e) such that 


p('=)(« 2 e 2 e FSie, K, m,V{e))) = ^ + Oe W, for k = {1,2} or {2}. 

P2 

However the presence of the first mutant speeds up the invasion of the second one: there exists a positive 
constant c such that for every e > 0; 

< (1 + «)!| 5 ) = f + 


^20 


20 


where S 20 > S 20 is defined by 


1 1 Ui ^ 1 

^20 S^\Sio “/V 5 / 52o’ 


where S = S '21 in case (2.71 and S = S '201 in case (2.8), and 

log if 


^ < <1+“>lr) = f - 7sf)+ 



Figure 2: Invasion times of the mutation 2 without and with the first mutation. The y-axis corresponds to 
the three type population sizes (0 in black, 1 in red and 2 in green), and the x-axis to the time. In these 
simulations, K = 1000, (/3i,(5i) = (2,0),i S S, Co,o = 1-8, Cop = 4,Co.2 = 3,Ci,o = l,Ci,i = 2.3,Ci.2 = 
3, C 2 ,o = 1.5, C 2 ,i = 1, C 2,2 = 2.1, and a = 0.5. 


2.3 How does clonal interference modify the invasion probability? 


Previous models predicted that clonal interference could only decrease the invasion probability of a mutant 
El- It was a direct consequence of the fact that the fate of a competition between two mutants was 
only dependent on the relative values of their fitnesses. Moreover, this implied that individuals with 
different phenotypes could not coexist for a long time. In our model, both decrease and increase of 
the invasion probability due to clonal interference may occur, depending on the competitive interactions 
between individuals and the time of appearance of the second mutation. Moreover, long term coexistence 
of several beneficial mutations are allowed when we do not assume that the mutant fitnesses are totally 


ordered. Recall the definition of invasion in (1.3). Then we have the following result: 


Proposition 3. Let Assumption{^ and one of the following conditions hold: 


S'oi > 0, 
S'oi < 0, 
S'oi < 0, 


S 201 < 0 and Assumption^ holds, 

S 21 < 0 and Assumption^holds 

S 21 < 0, S 20 > 0 and Assumption{^holds 
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Then for every n € \ { 0 } and (3,V > 0 ; 

lim =0, 

K^oo 


whereas there exist n G \ {0} and /3,V > 0 such that: 

lim p(2)(/(2,if,n,/3,y)) = §^. 

Notice that even if we recover here a classical result, saying that a mutant can get extinct because of 
the competition with an other mutant, we do not require 5'io > S '20 > 0: which would be the equivalent 
of assumptions done in population genetic models. 

Proposition 4. Let Assumptions^^ hold, and suppose that S 20 < 0. Then under one of the following 
additional conditions 


S'oi < 0, S '21 > 0 and < 


0 < a — ^ and < or 


„ - ^ > -A _ L. 

SlQ |SoiI S21 ’ 


5'i2 > 0 
r 

S 12 < 0 and S 02 < 0 


(2.9) 


Soi > 0 and S 201 > 0, (2.10) 

there exist n G K+ \ {0} and f3,V > 0 such that: 

K^oo 1^2 Pi 

where S = iS'21 in case (2.9 1 and S = S 201 in case ( [2l0t , whereas for every n G \ {0} and j3,V >0.' 


lim p(2)(/(2,/s:,n,/3,y)) = 0. 

if—foo 


Concerning the interplay of invasion and clonal interference, let us mention recent works which have 
taken into account the case where many benehcial mutations occur before any can hx [271[ll[Ti[2i[30]. 
The authors still assume transitivity of mutant fitnesses, but consider a regime of frequent mutations (high 
mutation rate or very large population). New mutations constantly occur in individuals already carrying 
other mutations, in their way of invasion, and the fate of a mutation depends on the genetic background 
of the individual where it occurs more than on its intrinsic advantage. They argue that this dynamical 
equilibrium is a way to preserve genetic diversity despite clonal interference, where the amount of variation 
results in a subtle balance between selection, which reduces it, and new mutations, which increase it. This 
approach is interesting and relates on experimental data which confirm that populations with so frequent 
mutations do exist |30j . However, due to the transitivity assumption, the authors need to assume that 
a large number of mutants co-occur in order to explain the possible coexistence of several types in the 
population, which is not necessary in our model. 


2.4 When beneficial mutations annihilate adaptation? 

An other interesting phenomenon can happen in our model: the occurrence of the second mutation can 
annihilate the effects of the first one and lead to the final Hxation of the wild type population 0 which 
would have been outcompeted by the first mutation alone. Such a phenomenon is also impossible in case 
of transitive fitnesses, as in this setting 5'io > 0 and «S' 2 i > 0 necessarily imply S 20 > 0. Proposition is 
illustrated in Figure 
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1 


Proposition 5. Let Assumptions^ and^ hold and suppose that 

So2 


Soi < 0 , Si2 < 0 < 5 ' 2 i , S 20 < 0 < 502 and 


|5i2||5oi| 5io ^21 |5oi| 


Then for every e > 0, there exist two finite constants /3(e) and V{e) such that 

p(i’^)(noeo G FS{e,K,/3{e),V{e))) = + 0^(e). 

Pi P 2 



Figure 3: Illustration of Proposition The y-axis corresponds to the three type populations sizes (0 in 
black, 1 in red and 2 in green), and the x-axis to the time. In these simulations, K = 1000, {Pi,Si) = 
(2,0),i G {0,1,2}, Co,o = 1-8, Ci^o = 1^2,1 = l,Ci^i = 2.3,Ci^2 = 2(70,1 = 5,(72,o = 2(7o,2 = 3,(72,2 = 2.1, 
and a = 1.9. 


2.5 At the origin of the Rock-Paper-Scissors cycles 

Rock-Paper-Scissors (RPS) is a children’s game where rock beats scissors, which beat paper, which in 
turn beats rock. Such competitive interactions between morphs or species in nature can lead to cyclical 
dynamics, and have been documented in various ecological systems [1I11I1I1111I1H1IS1I3S]. Let us describe 
two examples of such cycles. The first one |41j is concerned with pattern of sexual selection on some male 
lizards. Males is associated to their throat colours, which have three morphs. Type 0 individuals (orange 
throat) are monogamous and very aggressive. They control a large territory. Type 1 individuals (dark- 
blue throat) are polygamous and less efficient in defending their territory, which is smaller, having to 
split their efforts on several females. Finally type 2 individuals (prominent yellow stripes on the throat, 
similar to receptive females) do not engage in female-guarding behavior but roam around in search of 
sneaky matings. As a consequence of these different strategies, the type 0 outcompetes the type 1, which 
outcompetes the type 2, which in turn outcompetes the type 0. The second example [28] is concerned 
with the interactions between three strains of Escherichia coli bacteria. Type 0 individuals release toxic 
colicin and produce an immunity protein. Type 1 individuals produce the immunity protein only. Type 
2 individuals produce neither toxin nor immunity. Then type 0 is defeated by type 1 (because of the 
cost of toxic colicin production), which is defeated by type 2, (because of the cost of immunity protein 
production), which in turn is defeated by type 0 (not protected against toxic colicin). 

Neumann and Schuster |36] modeled such interactions by a three dimensional competitive Lotka- 
Volterra system. In particular they proved that migrations or recurrent mutations were not necessary 
ingredients to obtain limit cycles, as it was assumed in previous models (see nmsi for example). They 
studied the long time behaviour of the system but payed little attention to initial conditions, only assuming 
“the presence of all three strains in one homogeneous medium”. But the question of initial conditions is 
crucial. Indeed, how to explain the appearance of such cycles whereas when only two strains are present 
one of them is outcompeted by the other one and disappears? Our simple model provides a framework 
explaining how a cyclical RPS dynamics emerges in an ecological system thanks to the interplay of two 
successive mutations. Proposition [^ is illustrated in Figure]^ 
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Proposition 6. Let Assumptions^ and^ and the following inequalities hold: 

Sqi < 0 < Sio, Si2 < 0 < S21, S20 < 0 < S02, 

and 

( — 


0 < a — < min | 

Jio 


5, 


02 




02>J10 


I'S'oil’ |5 'i2||5'oi| ’ |5'i2||S'oi||S'2o|' 


1 

<5'21 


( 2 . 11 ) 


( 2 . 12 ) 


Then for every Z € N, i/ we eall “eycle” the interval between two local (non null) maxima of the type 1 
population between which type 2 and type 0 populations also hit one loeal (non null) maximum, and N the 
number of cyeles 

lim = 

K^oo Pi P2 

Moreover, if we denote by Di the duration of the kth cycle and introduce: 

1^01115: 


T{a,l,K) := (a- 7 ^ + 7 ;^Vl+ 


no 


s. 


21 


S, 


02 


Sn2S 


02'JIO 


n ^ ^ *^02*^21 * 5 in ^ 


then Dk satifies for a finite c 


(1 _ p K)<Dk<il + ce)T{a, l,K))=^^+ Ok^e). 


Pi P 2 



Figure 4: RPS cycles. The y-axis corresponds to the three type population sizes (0 in black, 1 in red 
and 2 in green), and the x-axis to the time. In this simulation, K = 1000, {Pi, Si) = {2,0),i G {0,1,2}, 
Co,o = Ci,i = 1 ^ 2,2 = 2, Co,i = 2.5, Co ,2 = Ci,o = C 2 ,i = 1, Ci ,2 = C 2 ,o = 3, and a = 1.1. 


Remark 2. IFe say that System (1.151 is permanent if there exists a compact attractor B C mt(]R{_) of 
its solutions, whose basin of attraction is mt(M{_). Under Condition (2.11), Theorem 2 in ]80( states that 
(1.15) is permanent if and only if 

|5'oi||5'i2||S'2o| < <S'o2>S'21*S'lO. 


This condition is satisfied under the assumptions of Proposition Hence if one of the types gets extinct, 
this is due to the demographic stochasticity and not to the behaviour of the approximating dynamical 
system. 

Moreover, we are not able to know in general if the interior fixed point is globally attracting of if there 
exist stable periodic orbits for the flows of the three dimensional deterministic Lotka- Volterra system. In 
Appendix we give two examples of systems satisfying the conditions of Proposition but with distinct 
long time behaviours. 
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3 First phase 

In Sections to we describe the dynamics of the successive phases. For sake of readability, we do not 
indicate anymore the superscript of the probabilities. The context will make clear the mutations which 
occur. The first phase is rigorously defined as follows: 

Phase 1 := {t > 0,N^{t) + N^{t) > 1, sup N^{t) < lsK\}. 

Under Assumption]^ the second mutation occurs during the third phase. This corresponds to the case 
already studied in [7] and we will recall the outcomes in this case in Section 3.3 Under Assumption or 
l^the second mutation occurs during the first phase and we have to study the resulting dynamics during 
the first phase. 


3.1 Assumption]^ 

Let i, j and k be distinct in {0,1,2}, e > 0 and iF £ N. We introduce a finite subset of N containing the 
equilibrium size of a monomorphic i-population. 




T^l - r, ^i,j + C*! 

Kirii — 2 £-— - - 


a. 




r, C'i J + C'i- 
^ 1,1 


nN, 


(3.1) 


and the stopping times Ta^’^\ /p{K,i) ^ which denote respectively the hitting time of size [aj 

for a £ M+ by the population of type i and by the total population of types i and j, and the exit time of 
jjy population of type i. 


:=inf{t>0,7Vf(t)= Laj}, 

:= inf [t > 0, N^{t) + Af (t) = [a\}, 
:= inf 


(3.2) 

(3.3) 

(3.4) 


Finally, we introduce a finite subset of N which may contain the type 2 population size at the end of 
the first phase. 




AT" 


?V‘S'2o(s}77-a+£) 






nN, 


(3.5) 


where we recall that a log K is the time of occurrence of the second mutation, and the invasion fitnesses 


have been defined in (1.6). Then we have the following possible states with positive probability at the end 
of the first phase. 

Lemma 3.1. Under Assumptions^ and^ there exists a positive constant Mi such that 

AK,i) / rj.{K,i) ^ f{K,0) y- p ql.'l ^ ^ 

/?! 

(6) ^ ^ = OK{e), I^je {1,2}, 

Pi ^ Pj ' 

(c) Af (T^^’'^) £ ) = If ^ + Oic(e). 


(a) P(ro^"’‘^ < W £ {1,2}) = (l - + OK{e), 


Proof. In the vein of Fournier and Meleard [20j we represent the population process in terms of Poisson 
measures. Let Qo(ds,d0), Qi{ds,d0) and Q2{ds,d0) be three independent Poisson random measures on 
with intensity dsd0, and recall that {ei,i £ {0,1,2}) is the canonical basis of Then the process 
can be written as follows: 

1 


A'^(t) = A'^(O) 


2=0 


fo 


le</ 3 iArf (s-) ~ (s-)<(if (Af^(s-)) 


+ 62 lt>a logic 


(-/ 

J CK 


Oi log K 


le</32Wf (s-) “ '^Q<e-P2N^ (s-)<d^ {N^ (s-)) 


Qi{ds,d0) 

Q2{ds,d0)^, (3.6) 
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where € {0,1,2}) have been defined in (1.2). The idea is to couple the population process with 

birth and death processes to get bounds on the different hitting times. To this aim, let us introduce 
approximations of the so called rescaled invasion fitnesses S^jfii, for i G (1, 2} and e > 0 small enough: 


s-o’ ^ := ^ {j3i - Ci^o + Se 


C'o,i + C[ 


0,2 


C( 


0,0 


) - {a,i + a,2)e) 


A 

These real numbers satisfy for s small enough 


(e.+) _ 1 /a (- o C'o,l + C' 0,2 {{ 

,0 5 , a.o(no 2£ jj. 


0 < s-o’ ^ ^ < 1, and 




’iO 


— S 


iO 


< 

“ A 


1 (4%^' 

^ c^o.o 


0,2 


+ Ajl + C'i.S )^' 


(3.7) 

(3.8) 

(3.9) 


Thanks to these definitions we can introduce, for * G {—, +} the processes 

7V(^’*)(t) = LnoKj + f I 
Jo J] 

^0<e-/3oJvA'*>(s-)<(5o+Co,oArA’*’('>“) + l{.=-}^(Co.i+Co,2))AfA’*’(s“)J^°^^®’'^^^’ 

and the supercritical birth and death processes 


_^e</3owA’*'(s-) 


= 1 + 


_le</3iivA’*’(s-) ^0<e-/3iArA'*^(s-)<(l-A''*>)/3iiv''’*’(s-) 


Qi(ds,d6»), (3.10) 


and 


>oc log K 


1 + 


> a log K «/M+ 




Then recalling Definition (3.2) we get 


Nt~\t) < N^{t) < a.s. V t < A 


Q 2 {ds, d9)^. 

(3.11) 

(3.12) 


and for i G (1, 2}, 

iVf’”^(t) < iVf (t) < Aff’+^(t) a.s. y t (3.13) 

Moreover by construction for every (*i, * 2 ) € { —, +}^ the processes and are independent. 


To prove Lemma 3.1 we first need to show that with a probability close to one, Couplings (3.12) and 


(3.13) hold during the whole first phase; to this aim we first prove the following asymptotical result: 

liminf P(A{^) > 1 — ce (3-14) 

iC —>00 


for a finite c and e small enough, where 


:= ^ r^iK,2) ^ j.{K,12) ^ fiK,0)y_ 


First applying (A.S) we get the existence of a positive V such that for * G {—, +} 

lim P(7;(^’*) > = 1, 

K—¥00 


where 

■= inf i 


(3.15) 


(3.16) 

(3.17) 
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We divide the probability in (|3.14|) into two parts according to the position of with respect to 

oVK. 


1 - P(Af) = A a t, 


(K,2) rp{K,12) 
eK '^-‘-0 


(3.18) 


+p(e^^ < A a 

Thanks to (3.12) and ( 3.16[ ) we get 

lim P(f;(^’°) < a a = 0. 

K—^oo 

Consider now the second probability in (3.18). The event < Te^'^Y means that Couplings (3.13) 

hold at least until time A A Hence 




a a C < (Tq^^’"^ V tY’~^) A A 


C {' 




<T, 


(A-) 


AT 


eif 


’}U{e 


<n 


A T 


( 2 .+) 


eK 


}• 


But thanks to Equations (A.2) and dwst we know that for i e {1,2}, 


lim P({e'^^ < TY’ A (Tg^*’ ■’= 00 })= lim Pde'^^ < tYk’} A {TY'^’ < 00 }) = 0, 

K—>oo K—>oo 

where A denotes the symmetric difference: for two sets B and C, B A C = {B O C^) U (C C B^- This 
implies that 


P(e'^^ < a Y 


{K,2) . rp{K,12). 


A 

eK ^^^0 


) < E A < 00 ) +0^(e). 


2 = 1 


But Definitions (3.10) and (3.11) imply that for i G {1,2} the event {Tq*’"'’^ < 00 ,T^' ^ = 00 } is empty. 
This ends the proof of (3.14). 


We are now able to prove Lemma 3.1 We assume that holds, which is true with a prob¬ 
ability close to one according to (3.14), and implies that Coupling (3.13) holds on the time interval 

[0,7AAi) 


A rp{K,2) . rp{K,12)-i 

'-eK '^-^eK -^0 J 


(a): Let us recall Definitions (3.7) and (3.8). Thanks to the independence of the processes and 

we get by applying (lA.l ) 


P(Tc 


( 1 .+) 


r^(l,-|-) rp(2, + ) 

^ ^eK Ao 


<T 


(2.-^) 


eK 


) = 






1 - 


1- (l-sl'^’+^)L^'^J 


1 - 


1- (l-s^‘^’+^)L^'^J 


> 


<^20 \ 


P 2 


for a finite c, K large enough and e small enough. Moreover we have the following inclusion: 

r aK rT,(l,-t-) ^ ^(1,-t-) T-i( 2, + ) ^ ry^(2,-t-)., r .K rp{K,l) rr>{K ,l) rji{K ,2) rji{K ,2) . 

l^e Ao ^ -‘^eK Ao ^ ^eK J l^e Ao ^ eK Ao < eK S' 

We then get 


P(T( 


(Ai) 


<T 


(iC,l) rr.(K,2) 


eK > P(Af,T, 


0 


^(1, + ) rp(2, + ) rTi{'^, + )\ 

^ -^eK 5-^0 ^ -^eK J 


'-eK 


> P(To^^’+) < < Tj^+)) - P((Af )d 


\ *^10 A *^20 \ 


S2 

A 


(3.19) 
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for a finite c, K large enough and e small enough, where we used (3.14| and (A.l). 


(b): The independence of the processes and again yields 


< T: 


,(!,-) ^(2,+) 
'^0 5^0 


( 2 .+) 

eK 


) = 






1 - 


1-V 1- 


\) - PA P 2 ) 


for a finite c, K large enough and e small enough. Moreover, thanks to Coupling (3.13) we get 
r ^ 7 ^( 2 , + ) , 7 ^( 2 ,+) 71 ( 2 ,+) ^ rp{l, — )'[ 

/ 4^ 7^(^4) ^ rp{K,l) rp{K,2) rp{K,2) rp(K,2) rp{K,l) ^ 

l^e Me_R- < 2g ,ig <■ ^eK Mo ^ ^ eK I' 


But according to Lemma A.l on the event {T^4 ^ ^ ^ finite constant 

( 2 .+) 


M such that for K large enough with a probability close to one < a log K + M, and ^ is close 

to log Ar/S' 4 ’~^- We finally get: 




A if 


> 

P(2A^.i) 

< 

^ 2g 

> 

p(j,(L-) 

< T^^’~ 

^ 2g 

> 

^(1- 

>520 \ 


Pi V 

a' ' 


^eK 

r(2.+ 

Cif 


^eK 

p(l.- 

Mif 


— ce 


(3.20) 


for a finite c, K large enough and e small enough, where we used (3.14) and (A.l). By interchanging the 
roles of 1 and 2 we derive a lower bound for ( 6 ). 


(c): Let us now focus on the last inequality. First using again independence between fV^’ ^ 3 'iid \ 
and (A.l I we get 

p(t4 


^eK 


'0 ^-^eK 


^( 2 .-)n . <S'io S20 

) > — - - -C£. 


° Pi P 2 

But as Coupling ( |3.13 l only holds before time A we have to determine which process, 

or hits [eA'J first. For i € {1,2}, from (A.2) and (A.3), 

P(AA < < 00) < C£. 


Using again (A.3) we get the existence of a finite constant c such that: 

>S'io S 20 


(rg(*’-)=oo,ie{i,2}) >^^ + 0;,(e) 


and 


•{{k 


=2}OL + 


1 — ce 


3i0 


log a: < U 


:i.+) 


n(i, —) 


eK ^ ^eK ^ 


(l{i= 2 }a! + “ 1 ^) log^:^o*’ ^ = 00 , ie { 1 , 2 }^ 




Hence, as under Assumption]^ l/S'io < 1 /-S' 2 o + Q;, processes (A'4’*\ * G {+, —}) hit [eATj before processes 
G {+j~}) on the event {Tg*’~^ = oo,i G {1,2}} with a probability close to one. The last step 
consists in determining the values of (A^4’*^ * ^ {+>“}) on the time interval [(1 — ce)log AT/S'io, (1 + 
ce)logA'/5'io]- First we notice that ( A.3|) implies for * G {+,—}: 


r 


( 2 .*) 


^320(^7^—2«) 


log a: 


’20 / 1 „ \ ni(2.*) 1 

-777Y U- a*2ce , a.s. on {Tg' '= 00 }. 

sTn’ ^-^10 ^ 


a (sf 
P 2 S 20 


C£ 
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Moreover, we get from (A.l) 


' rj^' 


( 2 -) 


k® 2 o(^- 


^-3ce) j^S2o{-i 


y.-ce) 




(1 - 


SonC-p-a —ce) / \ ^ 

V-*- *^20 ) 


-a-2ce) 


(i-4ro 




-( 1 - 




0, K ^ 00. 


’20 


This completes the proof of the lower bound for (c). Adding (|3.19|) and (|3.20|) ends the proof of Lemma 

KB □ 


We end the study of the first phase dynamics under Assumptions and by an approximation of the 
duration of this phase in the case we are interested in: 

Lemma 3.2. Recall the definition of Mi in Lemma \3.1\ Under Assumptions^^ an 


<( 


log a: 

1 - ce )—— < 
*^10 


< (1 + ce)^ iVf S ) > 1 + 0^(e). 

<^10 ' 


We do not detail the proof of this Lemma. It is a direct consequence of Couplings (3.13) and Equation 


(3.21). 


3.2 Assumption 

Under Assumptions and the type 2 population size has a positive probability to become larger than 
the type 1 population size during the first phase. Let us introduce a finite subset of N which may contain 
the type 1 population size at the end of the first phase. 


:= 


Sio Sio 


nN. 


(3.22) 


Then we have the following possible states at the end of the first phase. 

Lemma 3.3. Under Assumptions^ and[^ there exists a positive constant M 2 such that 


(a) 

(b) 

(c) 


^ ^ g 2}) = (1 - (1 - + OK(e) 


S 2 

P2 


OK{e), i^je {i,2} 




fil fi 2 


We do not prove this result as it is very similar to Lemma [3T| The idea is that when the 2-population 
survives the first phase, it takes a time of order log Ar/ 5'20 to hit the value [eATj, whereas the 1-population 
size needs a time of order \ogK/Sio to hit such a value. As under Assumption a -I- l/S '20 < l/5'io, the 
2-population size is the first to represent a positive fraction of the total population size. The value of the 
1-population size at the end of the first phase is obtained thanks to Coupling ( |3.13[ ) and Equation (A.3). 
We have also an equivalent of Lemma |3.2[ 


Lemma 3.4. Recall the definition of M 2 in Lemma \3.3\ Under Assumptions^ an d0 


“^(l—ce) ^ad- 


log AT 
S- 


20 


<T 


(K, 2 ) 


eK 


< (l+c£)(a+i^)|ri^’") < e = l+OKie). 


3.3 Assumption 

This case has already been studied in [7] and we said a few words about it just after the definition of 
Assumption!^ We here present rigorously the results of [7] with our notations: 
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Lemma 3.5. Under Assumptions^ anrf[^ 


^ J,(K,0)^ = I'l _ + OK(e), < Ti^’O)) = ^ + 


5i 


and 


1 — ce) 


Pi 


log if ^ ^(K,i) ^ , „^logif 


S 


10 


'^eK 


< (1 + ce)- 


s 


10 


rpi^d) - 
^eK ^ 




Pi 

= 1 + O/f (e). 


4 Phases 2 and 2n 


In this section, we describe the dynamics of ’’deterministic phases”, when some of the population sizes are 
well approximated by the solution of a two- or three-dimensional competitive Lotka-Volterra system. 


4.1 Two-dimensional case 

Let us denote by the end of the first phase, when at least one of the mutant population survives: 

rp(K,l) ^ rp(K,2) 


"eiC 


^eK ’ 


by i{Pi) the label of the first mutant population which hits the value [eifj, and by j{pi) the label of 
the other mutant population. We will focus on the most interesting case, when the population j{pi) does 
not get extinct, as the other ones have already been studied in [7] and [9], and introduce the event: 


B, 


K 


_ ^ )) rjl[K,0) AjK f rp{fe ))\ t( ^ J (‘f’f’) ) \ 

■— \^eK l(<^f) J ^ MiS j' 


(4.1) 


We will now consider the second phase of the sweep. It corresponds to the interval between the time 

when the mutant population i{p^) hits the value [eifj and the time when the rescaled pop¬ 
ulation process {N^, is close enough to the stable equilibrium of the dynamical system (1.7) 

with labels 0 and i{Pi')- To define rigorously the duration of the second phase we need to introduce a 
deterministic time tg {i,j) (see (4.3)) after which the solution of the dynamical system (1.7) with initial 
condition z = {zi,Zj) G is close to the stable equilibrium. To do that in a simple way we introduce a 
notation for the stable equilibrium of (1.7) independent of conditions (1.8) and (1.13): 


(eg) (eg) 

■ = 


Hj Cj , 
^1-3 




if ( 1 . 8 ) holds, 


n)7ej -I- niycj, if ( |1.13[ ) holds. 


where we recall Definition (1.14). Hence for e > 0, fe (hj) can be defined by 


:= inf{s > 0,Vt > s. 


|nf^(t)ei -f 


{t)ej - 


,(e9)| 




(4.2) 


(4.3) 


where ||.|| denotes the Li-norm on K^. As we do not know precisely the initial value of the rescaled 
process (N^, at the beginning of the second phase, we consider the supremum over the possible 


'■= {(^o, 2 :q 0 f )),\zo - no\ < 3e(C'o,i -k C'o, 2 )/C'o.o, e/2 < z^{pi ) < e}. 


and 


(4.4) 


f,(0, *(</f)) := sup{4^)(0, i(/>f)), z e }, 

which is finite. We are now able to define rigorously the second phase: 

Phase 2 := < t < -k 4(0, f(^f))}. 

During the second phase two populations have a size of order K, the wild type population and the 
mutant population of type {{pp). We will prove that the dynamics of these two population sizes are well 
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approximated by the deterministic two-dimensional competitive Lotka-Volterra system (1.7), which stays 
in a neighbourhood of its stable equilibrium after the time ^^(O, The duration of the second phase, 

te(0, )), does not tend to infinity with the carrying capacity K, unlike the durations of the first and 

third phases. As a consequence, the size of the )-population stays negligible with respect to K during 
the second phase. Recall Definition (4.1). Then we have the following result: 


Lemma 4.1. Recall the definitions of the Mi and M 2 in Lemmas \3.1\ and \3.^ respectively. Then for 
every e > 0, 


lim +s) e ^ 2 ^ < te(0,i(^f)), 

K—voo \ 


— Co -I- + ts{Qfi{4>i ))) “ 


'Oi(0f) 


<£2 A inf 


zeBi^ 


= 1 . 


Remark 3. Here we cannot apply directly (A.5) because it requires positive lower bounds for the jump 
rates of each population divided by K. Indeed in our case, the jump rate of the population 2 is of order 


which is negligible with respect to K. Hence we will couple the process and the process {N^, 

with well known processes to get bounds on their dynamics. 

Proof of Lemma 4-1 For sake of simplicity we will write {i,j) instead of {i{(j)^), j{(j>^)) along the proof. 
Let us first prove that during the second phase the process does not evolve a lot. First, notice that 
the processes Ni^ and are bigger than if they were evolving alone. Hence if we introduce the event 


Cf:={ sup {N,^{s)+N,^{s)}>2{no + n,)K}, 

^ s<te{ 0 ,i) ^ 


(4.5) 


we deduce from (A.4) that, 


lim 

K—^oo 


sup J 


) = 0 


(4.6) 


where we used the convention 


F(m„n,){.) ■■= P(.|(lVf, Af )(0) = (A„ A,)), (*,j) e fA 

To control the number of j-individuals during the second phase, we introduce two birth and death processes, 
A-*^’^’ ^ and constructed with the same Poisson random measure Qj as Nf- (see representation 

g§) 

)y(A2,-)/,N — _ / / i O (dsdO') 




0 JR+ 




Recall Definitions (4.1) and (4.5). We see that on the event D Cf, 

ArjA 2 .-)(^) ^ (<^f +t)< Aj^’"’+^(t), a.s. Vt < D(0,f). 

Let us first focus on the process It is a pure death process with individual death rate: 

S := 6j 2(ho -l- fii)(Cjfi Cj,i), 

and we can construct the following martingale associated with this process: 


Mf{t) := A^’^’”)(t)e^‘ - A^’^’”^(0) = - 


Qj{ds,d6), 
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where Qj is the compensated Poisson measure Qj{ds,d0) := Qj{ds,d9) — dsd9. Moreover, its quadratic 
variation can be expressed as 


{M^)t = 




^dsde= f 6N^^’^’~\s)e^^^ds= [ 5Mf{s)e^^ds. 

Jo Jo 

Then Markov Inequality leads to 


E 


< 


(^f)te(0,d 




_ e-^‘e(0.i)^ 

which implies that 

lim Pf inf {nI^’'^’~\s)} < 

K^oo Vs<tj(0,j) ^ ^ 

as a death process is non increasing. In the same way we prove that 

lim pf sup =0. 

K^oa \ s<t,(0,i) ^ 


^gitelO.z) _ 

^I _ -) 


2 > 


= 0 , 


From (4.6 1 to (4.71 we deduce that: 

lim pfivf (<^f + s) G < t,(0, *) ) 

K^oo \ ^ / 


= 1 . 


(4.7) 


(4.8) 


Now we want to control the dynamics of populations 0 and i during the second phase. We introduce 
two pairs of processes, and whose dynamics are well know and 

such that with high probability, 

< ivf < fcG{0,*}, 

during the second phase. These processes are defined as follows for t > 0 and (a, b) = {i, 0) or (0, i): 

' iVi^’"’-)(t) := iVf (0f) +/o7e+ 


iVf’2’+)(t) := TVf (</)f) + Q,(ds, dO) 




■'■0<e-/367v)^-"'+>(s-)<(56+C6,aAfl^'"--'(s-)+C6,tJv)^'"'+>(s-))Ar)^-^-+>(s-) 

Notice that conditionally on the initial condition of the second phase, the processes k G 

{0, i}, * G {—, +}) are independent of the process , and that on the event 


< oo} n {Nfi(l>i + s) G J 2 m^s^'Js < te(0,f)}. 


we have 


\s) < Nj^{(j)i + s) < a.s. Vt < te(0, j) and fc G {0, i}. 

Moreover, a direct application of ( |A.4 ) leads to 

lim sup P( sup - {n[^\t),n^^\t))\\> 5) = Q, 

(zo.Zi)64^’“V*^x[e/2.e] 0<t<ie(0.i) 

for £ and 5 > 0, where 5 denotes the complement of * in {—,+}, has been defined in (O and 

te(0, i) in (4.4). Adding (4.8) completes the proof of Lemma [4.1[ □ 
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The dynamics of the population process is a succession of phases where at least two types of populations 
have a size of order K (2nth phases, n G N) and of phases were at most one population type has a size of 
order K ((2n + l)th phases, n G Z+). The following lemma completes the description of the dynamics of 
phases 2n, in the case where two populations have a size of order K. It is a generalization of Lemma [dJ] 


and we do not give the proof. Recall Definition (4.4|. Then 


Lemma 4.2. Let {i, j, k} = {0,1, 2}, s > 0, n,Ci >0, 7 S (0,1) and assume that 


{N^ {0),Nj (0),7Vfc {0)) = with \n - n^\ < CiS. 


Then 


lim 

K—¥00 


(s) e < Ui.j). 


— {Ni + ej)%(i,j))-n,j 


(e?) 


< A inf 
zGBA 


n[''\te{i,j)) 


= 1 . 


4.2 Three-dimensional case 

From a probabilistic point of view, the case where the three types of populations have a size of order K 
is simpler, as we can approximate the rescaled population process by the three-dimensional deterministic 
Lotka-Volterra system (1.151 according to Equation (A.41. But the behaviour of the solutions of the three- 
dimensional Lotka-Volterra systems are much more various than these of the two-dimensional systems. 
They have been studied in detail by Zeeman and coauthors [48l|47l|46] and we will now present some of 
their findings. 

As in the case of two dimensional systems, the invasion fitnesses 


{Sij, Sijk, i,j, and k distinct in { 0 , 1 , 2 }) 

will determine the overall behaviour of the flows. In the two dimensional case for interacting populations 
of types i and j, there are three possibilities: 

• Either Sij < 0 < Sjf, then the only stable fixed point is njCj 

• Or Sij, Sji > 0; then the only stable fixed point is n^ij &i + 'hij'^^j 

• Or Sij, Sji < 0; then there are two stable fixed points, UiCi and njCj. 

In the three-dimensional case, Zeeman [48] numbered 33 equivalence classes. An “equivalence class” 
is a given combination of invasion fitnesses signs modulo permutation of the indices. Flows of systems 
belonging to the same class have the same type of long time behaviour, which is well determined in the 
25 first classes and more complex in the 8 remaining classes. These behaviours are represented in Figure 
and we will now explain their meaning. By an application of Hirsch’s Theorem |23j Zeeman proved 
that there exists an invariant hypersurface of denoted E, such that every non zero trajectory of a 
three-dimensional competitive Lotka-Volterra system is asymptotic to one in E for large times. E is called 
the carrying simplex, beeing a balance between the growth of small populations and the competition of 
large populations. A locally attracting fixed point is represented by a close dot •, a locally repelling one 
by an open dot o, and a saddle by the intersection of its hyperbolic manifolds. 

The classes 1 to 25 have no interior fixed points, and their dynamics are well known: the system 
converges to one of the stable hxed points with at most two positive coordinates. The solutions of systems 
in class 32 converge to one of the monomorphic equilibrium density hi for i G {0,1,2}. The solutions 
of systems in class 33 converge to the unique interior fixed points (see |17| for the two last assertions). 
In these classes 26 to 31 the system can have periodic orbits depending on the value of the parameters 
(see Figure]^. More precisely Zeeman proved that there exists systems with and without periodic orbits 
in each of these classes. We have no general criteria to discriminate between cyclical and converging (or 
diverging) behaviours of the flows in these classes. Note however that Hofbauer and Sigmund (Theorem 
15.3.1 in |24p and Zeeman and Zeeman (Theorem 6.7 in |46p provided sufficient conditions to have a 
global attractor (or repellor). 
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Figure 5: The phase portraits on E. A fixed point is represented by a close dot • if it attracts on E; by 
an open dot o if it repels on E, and by the intersection of its hyperbolic manifolds if it is a saddle on E. 
(figures found in [35] and modified to include results of |3Z] about classes 32 and 33) 



Figure 6 : Phase portraits in classes 26 to 31 in absence of periodic orbits, (figures found in [35j) 


5 Phases 3 and 2n + 1 


At the beginning of the third phase there are two possible cases: either there are already individuals of 
type 2 in the population; this corresponds to Assumptions or Then two types of individuals have 
sizes close to (defined in ( |4.2| )), and the last type has a size of order K'^ for some 7 S (0,1). Or 

the second mutation has not occurred yet (Assumption and there are only individuals of type 0 and 
1. Let us first focus on Assumptions and [^ Then at the beginning of the third phase we have, with a 
probability close to >S'ioS' 2 o/(/ 3 i/? 2 ), two kinds of initial conditions: 

• Either < 0 < then is close to equals \nK\ for some 

n e [inf,,gBif n[,^^(<£( 0 ,i(()if ))),sup^gB|f and belongs to 


Or > 0; then is close to AT,AT) and belongs to 


J2M 


K~.S- 


) 

In fact such initial conditions will also be found in phases 2n + 1, with n >2. Indeed we will see that 
in some cases there will still be two populations with sizes of order AT and one population with a size 
of a smaller order after the two first alternations of stochastic and deterministic phases (see Figure]^ for 
instance). Lemma 5.1 describes the dynamics of a stochastic phase with such initial conditions. Before 
stating this lemma, we introduce a finite subset of and a stopping time: 










) - (%■ e, +n)f ej) 


< 2£: 


\CijCjk — CikCjj\ + \CjiCik — CjkCiil 


\CiiCjj CijCj^i 


22 























































f{K,^3) > 0, (iVf (t), N^{t)) i 


Lemma 5.1. Let us take i, j and k distinct in {0,1,2}. Assume that Sij > 0 and N^{0) = \ for 

some 0 < 7 < 1 . 

• df Sji < 0, Ski > 0, A}*'(0) G [{fii — e‘^)K, (n^ + e'^)K] and N^(0) = ISK\ for some 0 < 6 < e"^, 

^ rj.{K,k) ^ = i^OK{.e), fo- ^ ^ 1 ^ 




^ Sk^ ’ 

(l-T-e)|S,il 


Sfci 


) = l+Oif(e), for 


. If 5,, > 0, Sk^J > 0, iV,^(0) G [{n^ - e^)K, (n^] + e^)K] for I G (i, jj, 

p(j 7 K,/c) ^ = i + OKis), 

• If Sji < 0, Ski < 0, (0) G [(n^ — e^)K, {ui + e^)K] and N^(0) = ISK\ for some 0 < 6 < e"^, 

^ rj.(K,k) ^ f{K,i)'\^ = i + OK{e), 

• If Sji > 0, Skij < 0, iV;^(0) G - e'^)K, +s‘^)K] for I G {i,j}, 

p^T^K,k) ^ f(K,^3)'^ =i^OKis). 

The proof of this Lemma is very similar to the proof of Lemma |3.1| It consists in comparing the 
population process with birth and death processes, which can be subcritical in this case, and to show that 
the equilibrium sizes [hiiGj or are not modified much by population(s) with size(s) 

smaller than [eil'J. Hence we do not detail the proof. We are also able to get approximation of the 
stochastic phase duration: 

Lemma 5.2. Let us take i, j and k distinct in {0,1,2}. Assume that Sij > 0 and N^{0) = \ for 

some 0 < 7 < 1. Then there exists a finite constant c such that, 

• If Sji < 0, Ski > 0, Nf" {0) G [(hi — e^)K, (hi + s'^)K] and N^(0) = ISK\ for some 0 < 6 < 

Tsk'"^ < = 1 + Oic(£). 


1 1 — 7 

> 


Ski 


1 - ce) log K < < (1 + ce) log K 

^ki ^ki 


. IfS,, > 0, Sk^J > 0, fVi^(O) G [in‘^ - e^)K, (n‘ff+s^)K], I G {r,j}, 


P((l - ce)^logK < < (1 + ce)^ \ogK = 1 + OK{e). 

^ ^kij ^kij ' 

If Sji < 0, Ski < 0, Nf^iO) G [(hi — e^)K, (hi + s^)K] and N^(0) = 16K\ for some 0 < 6 < e"^, 

p((i-c£)(^v^) logiP < < (i+ce)(^v^) 


= 1 + Ok{£)- 


If Sj, > 0, Sk^J < 0, fVi^(O) G [{n^ - e^)K, (h^^ + e^)K] for I G 

T{K,k) ^ ^ ^ 


1 - ce) i-^ log K < < (1 + ce) log K 

\^kij\ \^kij\ 
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Let us now describe the dynamics of the third phase when Assumption holds. As in Lemma 5.1 


there are five possibilities, depending on the signs of the invasion fitnesses and of the value of a. To take 
into account the state of the population at the end of the first phase, we introduce the following notation 
under Assumption]^ 


:=P(.|W,iVf )( 0 ) = (iVo,fVi) 


and the second mutant appears at time {a — 1/S'io) log AT) 


Recall that on the event < oo}. Lemma 3.5 ensures that (pf- is of order log A'/S'io with a probability 
close to 1. Then we can state the following properties for the population dynamics during the third phase: 

Lemma 5.3. Under Assumptions^ and^ 

• If Soi < 0, S '21 > 0 and l/|S'oi| < a - l/^io + l/5'2i, 

inf + OKis), 

• If ^oi < 0, S '21 > 0 and l/|S'oi| > a - l/^io + l/5'2i, 


inf 


Vmno,n,) < T^k"'’ a 


^l-|Soi|(a-s^ + ^+£)^^l-|Soi|(a-^ + ^-£) \ ^ ^ 

1 I P 2 


If Sqi > 0 , S'201 > 0 , 


inf 




. If Soi < 0, .521 > 0, 


No<e^K,NieI, 

• If <5oi > 0 , S '201 > 0 , 
inf 


inf Va<oo, 

r ^ ^ P 2 


||(lVo.Wi)/iC-(HS'.4i’)ll<^ 




= 1 - 


3201 


+ Ok{,s), Vo < 00 . 


Once again the proof is very similar to the proof of Lemma |3.1[ and we can also derive approximations 
for the total duration of the three third phases: 

Lemma 5.4. Under Assumptions^^ ond]^ there exists a finite constant c sueh that, 

• If Soi < 0 , S 21 > 0 and l/|Soi| < a — 1/Sio + I/S 21 , 


P((l-ce)( —+ ^) logX < < (l+ce)( —+—) logi^lTo^""’") < = l+Oi^(e). 

If Soi < 0 , S21 > 0 and l/|Soi| > a — 1 /Sio + I/S21, 

P((l-ce)(a+f-) logX < < (l+c£)(a+f-) logif|Ti^’"i < = l+OK{e). 


• If Soi > 0 , S 201 > 0 , 

P((l - ce)(a + f—) log a: < < (1 + c£)(a+ f—) log AT 

^ ^ *^201 ' ^ *-?9ni / 


9201 ^ 


rp{K,2) 

^eK ^ 


T™) =l + OK{e). 
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6 Proofs of Propositions 1 to 6 

We now prove the main results of this paper. The proofs are direct consequences of Lemmas stated in 
Sections ii and[^ As the proofs of Propositions [l] and are very similar, we only prove the second one. 

Proof of Proposition^^ Let us first consider that Assumptions andand Condition \2.7\ are satisfied. 
(1 — 5'io//3i) 5 * 207 / 32 : According to Lemma 


2 -population survives the first phase and hi 


3.1 

isT 


with a probability close to (1 — S'io// 3 i)<S' 2 o //32 only the 
le size [eATj. Then it ’deterministically’ competes with the 


0 -population and outcompetes it as S '02 < 0. In this case the duration of the invasion, logAr/S' 20 , is not 
modified by the presence of the mutant 1. 


5 'io<S' 2 o/(/ 3 i/ 32 ): According to Lemma 3.1 and 3.2 with a probability close to S 10 S 20 /both 1- and 


2 -populations survive the first phase, which has a duration close to 

log Lf 


no 


The 1-population size is the first to hit \eK \, whereas the 2-population size belongs to (defined in 

(|3.5|)) at the end of the first phase. Lemma 4.1 and Markov Property imply that with a probability close 


to one, the 2 -population size stays of the same order during the second phase, and the 0 - and 1 -population 
sizes become close to {aK,niK), with 0 < o < depending on ecological parameters but not on K, at 
the end of the second phase. As S '21 > 0 , the 2 -population size grows during the third phase and takes a 
time close to 

to hit the size [eK\ (Lemma |5.2[ ). Furthermore at the end of the third phase the 0-population has a size 
negligible with respect to K (which can be 0, see Lemma [aT|) . During the fourth phase the 2-population 


’deterministically’ outcompetes the 1-population, as 5 'i 2 < 0 < 5 ' 2 i (Lemma 4.2 1 . Then the 0- and 1- 
populations get extinct, as S 02 and 5 'i 2 < 0 (Lemma |5.1| ). Such a trajectory is illustrated in the second 
simulation of Figure In this case the total duration of the mutant invasion is 


/ 1 

1 

f ^ / 1 





) - a) 


) logi^, 


(the —a log AT is due to the time of the second mutation occurrence) which has to be compared with the 
duration of the invasion in the absence of the first mutation, logAr/S' 2 o- The resolution of a quadratic 
equation leads to the condition S 20 s]S'io/(l — aS'io), >S' 2 i [. More precisely, the inequality 


/1 

1 

^ / 1 





) -«) 


log a: 


520 


is equivalent to 


— 

*^20 c 
*2 


(1 \ 

^ f 1 

1 N 

( 7.-O! ] 

— * 5^20 ( — 

H— ^ — — Oi- ) 

\Sio 7 


^21 7 


whose solutions are 5*20 = 5 ' 2 i and S 20 = >S'io/(l — aS'io). As under Assumptionsand 

‘S'lo 


1 1 

~c ^ ~c ^ 
<J20 <310 


<S'20 < 


1 — aS 


10 


this ends the proof of this case. 


1 — S 20 /fi 2 - Finally with a probability close to 1 — S 20 /1^2 the 2 -population gets extinct during the first 
phase (Lemma |3.1[ ). 

The proof of the second case follows the same ideas. The only difference is that when the 1-population 
survives the first phase, the 0- and 1-populations coexist during the third phase, as 5'io and Sgi > 0. This 
ends the proof of Proposition □ 
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Again the proofs of Propositions and are very similar, and we only prove the second one. 

Proof of Proposition^^ Let us first suppose that Assumptions and hold, and that £'20 < 0, 'S'oi > 0, 
*5201 > 0 . 


1 ~ 5io//3i: According to Lemma 3.5 with a probability close to 1 — Siq/I3i the mutants 1 do not survive 


the first phase. Then following (ATS) we get that the 0-population size stays close to its equilibrium uqK 


during a time of order where P is a positive constant independent of K. Hence it is still close to 
this value when the second mutant appears (time a log AT). The latter one has a negligible probability to 


survive as the invasion fitness S20 is negative. After the 2-population extinction, using again (A.51, we get 
that the 0-population size stays close to its equilibrium value UqK during a time larger than where 
P is a positive constant independent of K. 


*5io/,5i(l — 5201//32): According to Lemma 3.5 with a probability close to Sio/j3i the mutants 1 survive the 


first phase, and the 0- and 1-population sizes get close to their coexisting equilibrium ([ho^^ATj, L^oi^AiJ) 
during the second phase. Then we get from (A.51 that the 0- and 1-population sizes stay close to this 
coexisting equilibrium during a time larger than for large K where P is a positive constant indepen¬ 
dent of K. Hence they are still close to this state when the second mutant appears (time a log AT). The 


latter one has a probability close to 1 — S201/P2 to get extinct before hitting [eATj (Lemma 5.3). After 

the 2-population extinction, the 0- and 1-population sizes stay close to their equilibrium value 

and [hg^^ATj during a time larger than where P is a positive constant independent of AT (Equation 

(M). 


<S'io*S' 2 oi/(/^i/^ 2): Applying again Lemma 3.5 and Equation (A.5), we get that with a probability close to 


Sw/Pi the 0- and 1-population sizes get close to K, K) and are still in this configuration when 

the second mutant occurs. The latter one has a probability close to 8201/^2 to hit a size of order AT 


(Lemma 5.3) whereas (A^o, A^i) still belongs to Then the final state depends on the signs of the 


other invasion fitnesses. The approximating deterministic Lotka-Volterra system after the third phase can 
belong to 7-12, 29, 31 or 33ird class described by Zeeman (see Figures]^ andcase p^ ). In all the cases 
the density of 2-type individuals does not tend to 0 during a time larger than where P is a positive 
constant independent of K. 

The proof of the second case follows the same ideas. We just have to take into account the value of 
a log K compared to some invasion time to know whether the 0-population gets extinct before the hitting 
of \eK\ by the 2-population. We detail this kind of comparison in the proof of Propositionand end here 
the proof of Proposition □ 

Let us conclude this section with the proof of Proposition The proof of Proposition follows the 
same outline and is simpler; hence we leave it to the reader. 


Proof of Proposition^ Suppose that Assumptionsand Conditions (2.11) and (2.12) hold. 


1 ~ 810 /See the beginning of the proof of Proposition 
(<5'io//3i)(l - 821/I32)- Lemma 


3.5 


and Equation (A.4) imply that with a probability close to ^loZ/Si, the 


1-population survives the first phase and the population state at the end of the second phase satisfies 
{No,Ni){(j)^) = ([noATj, \ niK\) with 


inf 0 
zGBl^ 


n^''^(te(0, 1)) 


< no < and Ini — ni| < £, 


where 1^(0,1) has been defined in (4.4). Then with a probability close to 1 — 821^^2 the 2-population does 
not survive, and the 0-population size, which can be compared to a subcritical birth and death process 
as S'oi < 0, also hits 0 while the 1-population size is still close to [hiTLj (Lemma |5.3[ ). In this case, the 
1-population size stays close to its equilibrium value [hi77j during a time larger than where E is a 
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positive constant independent of K (Equation (A.5)). 


-S'io<S'2i/(/^i/ 32): First applying again Lemma 3.5 and Equation (A.4) we get that with a probability close 
to S'lo//?!, the 1-population survives the first phase and the population state at the end of the second 
phase satisfies {Nq, = ([noATj, [niK\) with 

• r 1)) 2 11 - I 

mf —2- < hq < and \ni — ni| < e. 

zGBf 2 


Then Lemma 5.3 implies that with a probability close to 821 / 1^21 the 2-population size hits the value \ sK\ 
before the extinction of the 0-population, because 


(s'ln liSm I) 


log K, 


"S'lo l-Soii 

which is the approximate extinction time of the 0-population (Lemma |5.4[), is bigger than 




log AT, 


which is the approximate hitting time of \sK\ by the 2-population size (again Lemma 5.4). Combining 
Lemmas |3.5[ |5.3| and |5.4| we even get that at the end of the third phase, the 0-population size is of order 


i^l“l'®oil(“-s5y+sk), 

and that the duration of the third phase is 

(a - log a:. 

\ Oio *321^ 

During the fourth phase we can approximate the dynamics of the 1- and 2-populations, which have a 


( 6 . 1 ) 


size of order AT, by the system (1.7) with {i,j) = (1,2) (Lemma |4.2| ). As S '12 < 0 < 5 ' 2 i, the 1- and 
2-population states at the end of the fourth phase are ([riiATj, [712 ATJ) with 

. 2)) 2 j I - I ^ 

mf --- < ui < and |n 2 — n 2 | < e. 

zGBf 2 

During the fifth phase, the 0-population size has an evolution comparable to this of a supercritical birth 
and death process (S '02 > 0) and takes a time 


|5oi 


5, 


02 


(a - log a: 

\ Oio O21 / 


( 6 . 2 ) 


to hit lsK\ (Lemmas 5.1 and 5.2). The 1-population size has an evolution comparable to this of a 
subcritical birth and death process (S '12 < 0) and takes a time 


log AT 

IS 12 I 


(6.3) 


to get extinct. As according to Condition ( |2.12 1 

1*5^/ 1 

?02 V 


“ Sio + S 


< 


1^121’ 


(6.2) and (6.3) imply that the 0-population size hits [eA'J before the extinction of the 1-population, the 
fifth phase has a duration of order 


|5oi|/ 1 

— [a —— 



log AT, 


(6.4) 
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and at the end of the fifth phase, the 1-population size is of order 


K 


1 - 


I‘5l2l l‘^01 I 

•502 


(cx— - \--q^) 

'' 5io ' S 21 '' 


We then again apply the same reasoning: during the sixth phase, the 0-population outcompetes the 
2-population. During the seventh phase, which has a duration of order 


1 ^ 12 ! |5o 


-(a - \ogK, 

5 ^ *^10 *^21 '' 


S 10 S 02 ^ Jio 021 - 

the 1-population size hits the value \eK\ and the 2-population sizes ends at a value of order 


(6.5) 


K 


|g20llSl2l|S0ll 1 _)_ 1 ) 

OinOm ' oin 001 ' 


where we used Condition (2.12) which implies that 

1^1211^0 


*S'io<5'o2 




1 1 \ 1 
“ 310 ^ 321)^^18201 


Adding the durations of the third, fifth and seventh phases in (6.1), (6.4) and (6.5), we get the duration 
of the first cycle, 

1 1 \ / . I^oil |*S'oi||5'i2| 


1-f 


log K. 


V Sio >S'2i/V 502 5 o25io 

Then by induction we prove that at the end of the phase 3 + I, I G N, the 0-population size is of order 

at the end of the phase 5 -I- /, Z S N, the 1-population size is of order 

1 ISi^lISm I f ,. 1 I 1 '|c lSiUISnil|S,,nl y 

Sq2 ^ 5io S 21 5io5o2521 ' ^ 

and at the end of the phase 7 -I- Z, Z e N, the 2-population size is of order 

1 _ |52ol|5i2l|5oil / _ 1 I 1 W |5i2 I I^qi 11^20 I \l 
7^ 5io5o 2 ^ 5io S21'''' 5 io5o2521 ^ 

We also prove that the 3 -I- Zth phase, Z G N, has a duration close to 

' |5i2||5oi||52oh ‘ 


3io3o232i 

the 5 -I- Zth phase, Z G N, has a duration close to 


/ \ OlO 021 / 


1 

21 ' 


^10^02^21 ^ ^02 
and the 7 + Ith phase, I G N, has a duration close to 


7 |5i2||5oi||52ohM5oil 7„ _ ^ ^ ^ 

( q q q j <? V S' S / 

^ 'J10>5o2>->21 ^ <302 ^ *510 >521^ 


/ |5i2||5oi||52o| y |5i2||5oi| / _ J_ , M , ^ 

V 5io5o252i J 3io3o2 V 5io ^217 ^ ' 


5io5o252i 7 5io5o2 
This completes the proof of Proposition 


□ 


A Technical results 

This section is dedicated to technical results needed in the proofs. We first recall some facts about birth 
and death processes. They are used in Sections to and can be found in [T] : 
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Lemma A.l. Let Z = {Zt)t>o be a birth and death process with individual birth and death rates b and 
d. For a S M+, Ta = inf{t > 0,Zt = [aj} and Pi (resp. Ei^ is the law (resp. expectation) of Z when 
Zq = i G N. Then 

• For {i,j,k) G such that j G (i,k), 


Pj(Tk <Ti) = 

• If 0 < d b, for every i G Z+ and t > 0, 


MTo <t)=[ 


1 - {d/by-^ 

1 - {d/by-^' 

d(l - 6(^-6)*) y 


(A.l) 


(A.2) 


• If 0 < d < b, on the non-extinction event of Z, which has a probability 1 — {d/b)^°, the following 
convergence holds: 


T^/logN -G {b — d) , a.s. 

N—¥tyo 


(A.3) 


We also need large deviation results to quantify the difference between the rescaled population process 
N^/K and the approximating deterministic Lotka-Volterra processes (1.7) and (1.15). The following 
statements can be found in [7] Theorem 3 (b) and (c) and in [3] Proposition A.2. They follow from Dupuis 
and Ellis [TB] (Theorem 10.2.6 in Chapter 10): 

Lemma A.2. Let C be a compact of (IRlj_)^ ^ {0} or (M^)^, and T a finite positive constant. Then for 
every positive 5, 


lim supPfsup||A'^(t)/A-n(^)(i)ll > ^ ^o^(O) = VzK\) = 0, 
K^oo zee ^t<T > 


(A.4) 


where [zAj = ([zqAJ, [ziAj, [Z 2 AJ), and is the solution of (1.15) with initial condition z. 

Let h denote a stable equilibrium of a competitive Lotka- Volterra system in dimension one, two or three, 
with all coordinates positive. Let £ > 0 and denote the population process with the same ecological 
parameters as the considered Lotka-Volterra system and carrying capacity K. Then there exists a positive 
constant V such that 


lim P( sup 

K^oo \ +<'^vi 


t<e^^ 


N^(t) 


K 


< 2 £ 


A-^(O) 


K 


< £ = 1 . 


(A.5) 


Let us now prove Lemma 2.1 which precises the conditions needed to have transitive interactions 
between several types of individuals: 


Proof of Lemma \2n\ Let i,j, k he in £ and recall notation (2.1). The relations 

i ^ j and j F k 

are equivalent to 


or in other terms 


—-Ci,-< 0 < ^-C,, and -Cjk <h <—-Ck], 
Pj Pi Pk Pj 


and —<CjkA^. 
Pk Ckj 


1. If (2.2) holds, then 



Po 



Pk( 

'fh _ 

C.fc) 

= Pk( 


'-Pk 

/ 

\ 



1 

\ ~ 

Pk( 

(c„ 

A ^ 

)[C,k 


%k) 


Ckj- 


-Y 

cj 
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and 


Ski — 


> 


J^-Ck?l=pJ^^-Ck?j 

\pi J \ Pj Pi / 

P^iiCkj V V ^) - Ck^) > Pfc((Ci V - C 2 ) > 0, 

^ij 


which implies that 


i <k. 


2. If (2.3) holds, then one of the two previous inequalities is not satisfied. Hence either i k or i = k. 


3. Let us now assume that (2.3) holds. We can choose p > 0 such that 




> Cl and 




< C 2 . 


Assume now that 


Pi Pj 

Pj Pk 




P- 


Then 


Pj 


Pj 


Pk 


Pi 


Pk 


—-C 2 <0< ^-Ci, ^-^2 <0< —-Cl and — - Ci > 0 > — - C 2 . 


Pj 


Pk 


Pj 


Pk 


This ends the proof of Lemma |2.1[ 


□ 


We now present two examples of three dimensional Lotka-Volterra systems satisfying conditions of 
Proposition 1^ and exhibiting different long time behaviours. In the first one, the solutions converge to a 
limit cycle with constant period whereas in the second one they converge to the unique globally attracting 
fixed point of the system. Before these two examples, we need to give a definition and recall a result of 
Hofbauer and Sigmund. 

Definition A.l. A matrix A = (aij)o<ij <2 G Al 3 (IR) is called Volterra-Lyapunov stable if there exist 
positive numbers do,di,d 2 such that 


Let us introduce the matrix 


Then we have: 


diaijUiUj <0, Vn ^ 0. 

0<i,j<2 

Co,o Co,i Co,2 

= I Ci,o Ci,i Cl,2 

C2,0 C2,l C2,2 


Theorem 1 (Theorem 15.3.1 in [24]). If —A is Volterra-Lyapunov stable then the Lotka-Volterra system 
(1.15) has one globally stable fixed point. 

This theorem is needed to construct the second example in the following illustration: 


Illustration of Remark^ 1. In |36j . the authors model the RPS interactions of escherichia coli strains 
by the following system: 

{ ho = {S - niriQ - pni - p,n2)no, 

hi = {p- pno - Kani - pn2)ni, 

«-2 = (/3 - (m + 7)^0 - pni - pn 2 )n 2 , 

where all the parameters belong to . They give an example for which the solutions converge to a 
limit cycle with constant period. It corresponds to the following values of the parameters: 

p = p = 1, Ki = P = 2, K 2 = 1, 75, 7 = 2,84, K 3 = 0,844.. and (5 = 1,156... 


30 









For these values of the parameters we have: 


S'lo = 0,422..., 5oi =-0,0288... 

S'21 = 0,816..., S'i2 =-0,0143 
5o2 =0,0131..., ^20 = -0,220.... 


These invasion fitnesses satisfy (2.11) and 

5 ' 2 i > I'S'oil) 5 'o 2 >S' 2 i > |S'oi||S'i2| and S'o 2 <S'io 5 ' 2 i > |<S'oi|| 5 'i 2 ||S' 2 ol- 

Hence we can choose a such that (2.12[) holds. 


2. Recalling the definition of invasion htnesses in (1.6) we get that —A is Volterra-Lyapunov stable if 
there exist positive constants d^^di, ^2 such that for every n ^ 0, 

~ (A -5^- > 0 , 


a 


Pi - 5 ^ Pj - 5 j 


0<iJ<2 

which amounts to the existence of positive constants c/q, d 2 such that for every n 7 ^ 0 , 

diiPi - Si - S^j)ninj > 0 . 

0<iJ<2 

Expanding the sum and using the notation pi = Pi — 6i for i G £ we get the condition 


doPo^Q + dipin\ + d2P2n\ + (io(po + I'S'oil) + di{pi — S'lo 


rfo(po ~ S02) + ^2(^2 + IS20I) 

Let us now make the following choice: 


no 712 + 


noni + 

di{pi + IS12I) + d2{p2 — S21) 


nin 2 > 0 . 


, , ■ ^ c Po — ^02 Pi — Sio P 2 — S21 

p,di = 1, 7 G S, - = - = - = 7j, 

Po Pi Pi 

where 0 < 77 < 1/2. The condition to satisfy becomes 


I'Soil IS12I IS20I 


Po 


Pi 


P2 


= V, (A.6) 


ng + ni+n^H—i;-^(2noni+2non2 + 2nin2) = 


1 “ 277^^2 , 2 , 2 i ,1 + 27?. , , .2 


-(no+ni+n 2 ) + - 


-( 770 + 771 + 772 ) > 0, 


2 2 • 2 

which holds for every non null n € . The end of the proof consists in noticing that we can indeed 

choose the parameters as in (A. 61 . Assume for example that the pPs and the Ci^s are given. Then 
it is enough to take 

C-i-i,* = 77 ^^C^i, Ci i+i = (1 + 77 )-^C'i+i_i+i, where 7 is modulo 2 . 

Pi Pi+l 

Finally it is easy to check that the conditions of Proposition are satisfied. Applying we get that 
the Lotka-Volterra system (1.15) has one globally stable fixed point. This ends the construction of 
the second example. 

□ 


B Complete description of possible dynamics 

We now give a complete description of the possible population dynamics. Sections |B.1| and m are 
dedicated to Assumptions and respectively. In Section |B.2| we explain how the dynamics under 
Assumptions and can be deduced from the dynamics of the process under Assumptions and 
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Case 

Duration of the sweep: log K. 

Final state 

[Aj 

1 

Sin 



sio + S20i(l “) 


S2i<O.Si2<0,So2<0, (0,0.fla) (8) 

S2i>O.Si2<0,So2<0, 012 (29) 

S2i>O.Si2<0,So2>0,Sio2>0 012 (31) 

S2i>O.Si2<0,So2>O.Sio2<0 (n^gf,0,n^gl'>) (10) 

S 21 >0,Si2>0,So2<0,Soi2<0 (OjniabnJa^) (9) 

Si2>0.So2<0,< '^ 2 i< 0 or 012 (29 or 31) 

y 'S’21>0,5o12>0 

S21>0,Sl2>0,So2>0, (ftoi2’"oi2>”oi2) (33) 


Sio S 21 Sio J 

(0.fi(y.n(|>) 

© 

sio + sii (1+ So'is ) ( sio “ 

) 

5102<0. ing\0,ng^) (10) 

Sio2>0. 012 (31) 

[E 



(0.0,n2) 


si^+sii (i+^s§|^) ( 

) 


© 

s' +s' ( 1 +'s‘’'' + 's°'s 

Sio S 21 OQ2 S 02 S 102 / \ 

3-a) 

10 / 

012 (29) 

[H] 


) 


0 

SI^+^ 


(0,0,fl2) 

0 

Sio S 21 Sio / 



_Sjn_ 

(O.ni.O) 


Table 1: Possible dynamics under Assumption Encircled letters correspond to these in Figure See 
Section B.l for the caption. 


B.l Assumption!^ 

In Figure and Table we assume that the 1- and 2-populations survive the first phase. This happens 
with a probability close to S'ioS'2o/(/3i/32) (Lemma |3.2[ ) and is the case we are interested in, as otherwise 
there is no clonal interference and we are brought back to the invasion of one mutant already studied 
in [7] and [3]. The behaviour of the population process after the first phase depends on the relations 
between the ecological parameters of the different individual types. We describe the different conditions 
which discriminate between different scenari in Figure and list them in their chronological order of 
appearance. We also indicate the phase during which the conditions have an impact on the population 
process behaviour. 

In Table we describe the ’’final state” of the population processes under the conditions [A] to [1^ . 
For sake of simplicity we call ’’final state” in this setting the element of 


FS:= p|{nGK3,3^(e),y(e) > 0,lim ini e FS{e, K, l3{e),V{e))) > 0}, 


£>0 


K—^oo 


(B.l) 


when this intersection is non empty, otherwise we call ’’final state” the long time behaviour of the three 
dimensional Lotka-Volterra system close to the rescaled population process once the three population 
types have a size of order K. In this case we indicate the corresponding class of the dynamical system in 
the Zeeman representation (Figure]^ and we write ”012” when the invasion fitness signs do not allow to 
discriminate between a cyclical or stable coexistence of the three types of populations. 

We have also indicated in Table [l] the durations of the sweep. They have to be understood in the 
sense of Lemma 3.2 for instance, when the final state FS described in (B.l I is non empty. They are good 
approximations of the durations with a probability one up to a constant times e. When FS is empty, 
these durations correspond to the time needed for the three populations to hit a size of order K. 
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Figure 7: Different dynamics under AssumptionSee Section fB.!! for the caption 

B.2 Assumption 

Let us now explain how we deduce the possible dynamics under Assumption from the possible dynamics 
under Assumption Applying Lemmas |3.1| to |3.4| we get the dynamics represented in Figure with a 
probability close to S'io5'2o/(/3i/32) under Assumptions (|^ and) and respectively. We deduce that to 
obtain an equivalent of Figure it is enough to: 

1. Interchange the roles of the 1- and 2-populations, 

2. Replace a by —a in the conditions of Figure and in the duration of the sweep 

3. Add a log K to the duration of the sweep. 

We see that we do not get new long time behaviours by comparison with Assumptions and This is 
why we did not detail this case. 


B.3 Assumption 

Figure]^ and Tableare the analogues of Figure and Table for Assumption]^ Here we only assume 
that the 1-population survives the first phase (probability close to 5'io//3i according to Lemma 3.11, as the 
second mutation occurs during the third phase. The captions are the same, except that we add in Table 
[^a final state ’’cycles Rock-Paper-Scissors” which corresponds to the class 27 in Zeeman’s classification 
(Proposition]^. 
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Figure 8: From Assumptionto Assumption(see Section B.2| 
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Figure 9: Different dynamics under Assumption]^ Same legend 


as in Figure]^ 
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Case 

Duration of the sweep: log K. 

Final state 


A 



("m -"oV .0) 




cx + — 

i>201 

5’i2 <0,S’o2 <0,5’2o >0, (0,0,fi2) (7 or 8) 

S2i>O.Si2>0,Soi2<0, (0,n[\\n[^^) (9 to 12) 

S2o>O.So2>0,Sio2<0, (9 to 12) 

Si2>O.So2<0,S2o>0, 012 (29 or 31) 

'S’ij>0,V(ij')e{0,l,2} , (tloi2’tloi2’ttoi2) (93) 


C 





A 


+ fa—+ 

•521 5oi2 Y 5io 521 J 

S20>0,Si02<0, (n(g°J,0,nQ^) (10) 

S2o>0,Sio2>0. 012 (31) 

S2o<0, 012 (29) 


E. 


52ni 

(0,0,n2) 


© 



S2o>0, (ri^jbo.ftJij') 

‘S'20<0, (no,0,0) 





("02’’0,"^’) 


'ii 



012 (29) 


0 


“+sk + '^(l+^) (““^+sk) 

cycles Rock-Paper-Scissors 




a+o^^— 

52m 

(0,0,n2) 


K 


a+o^^— 

52m 



[lJ 


_Sjo_ 

(0,ni ,0) 


Table 2: Possible dynamics under Assumption]^ Encircled letters correspond to these in Figure]^ See 
Section B.3 for the caption. 


36 




















References 


[1] K. B. Athreya and P. E. Ney. Branching processes. Dover Publications Inc., Mineola, NY, 2004. 
Reprint of the 1972 original [Springer, New York; MR0373040]. 

[2] N. H. Barton. Linkage and the limits to natural selection. Geneties, 140(2):821-841, 1995. 

[3] J. P. Bollback and J. P. Huelsenbeck. Clonal interference is alleviated by high mutation rates in large 
populations. Molecular biology and evolution, 24(6):1397-1406, 2007. 

[4] L. Buss and J. Jackson. Competitive networks: nontransitive competitive relationships in cryptic 
coral reef environments. American Naturalist, page 223-234, 1979. 

[5] D. D. Cameron, A. White, and J. Antonovics. Parasite-grass-forb interactions and rock-paper-scissor 
dynamics: predicting the effects of the parasitic plant rhinanthus minor on host plant communities. 
Journal of Ecology, 97(6):1311-1319, 2009. 

[6] P. R. Campos, C. Adami, and C. O. Wilke. Modelling stochastic clonal interference. In Modelling in 
Molecular Biology, page 21-38. Springer, 2004. 

[7] N. Champagnat. A microscopic interpretation for adaptive dynamics trait substitution sequence 
models. Stochastic Processes and their Applications, 116(8):1127-1160, 2006. 

[8] N. Champagnat, P.-E. Jabin, and S. Meleard. Adaptation in a stochastic multi-resources chemostat 
model. Journal de Mathematiques Pures et Appliquees, 101(6):755-788, 2014. 

[9] N. Champagnat and S. Meleard. Polymorphic evolution sequence and evolutionary branching. Prob¬ 
ability Theory and Related Fields, 151(l-2):45-94, 2011. 

[10] P. Collet, S. Meleard, and J. A. Metz. A rigorous model study of the adaptive dynamics of mendelian 
diploids. Journal of Mathematical Biology, page 1-39, 2011. 

[11] C. Coron. Slow-fast stochastic diffusion dynamics and quasi-stationary distributions for diploid pop¬ 
ulations. Journal of Mathematical Biology, page 1-32, 2015. 

[12] C. Coron et al. Stochastic modeling of density-dependent diploid populations and the extinction 
vortex. Advances in Applied Probability, 46(2):446-477, 2014. 

[13] T. L. Czaran, R. F. Hoekstra, and L. Pagie. Chemical warfare between microbes promotes biodiversity. 
Proceedings of the National Academy of Sciences, 99(2):786-790, 2002. 

[14] J. A. G. de Visser and D. E. Rozen. Clonal interference and the periodic selection of new beneficial 
mutations in escherichia coli. Genetics, 172(4):2093-2100, 2006. 

[15] M. M. Desai and D. S. Fisher. Beneficial mutation-selection balance and the effect of linkage on 
positive selection. Genetics, 176(3):1759-1798, 2007. 

[16] P. Dupuis and R. S. Ellis. A weak convergence approach to the theory of large deviations. 1997, 1997. 

[17] S. Ethier and T. Kurtz. Markov processes: Characterization and convergence, 1986, 1986. 

[18] R. A. Fisher. On the dominance ratio. Proceedings of the royal society of Edinburgh, 42:321-341, 
1922. 

[19] R. A. Fisher. The evolution of dominance. Biological reviews, 6(4):345-368, 1931. 

[20] N. Fournier and S. Meleard. A microscopic probabilistic description of a locally regulated population 
and macroscopic approximations. The Annals of Applied Probability, 14(4):1880-1919, 2004. 

[21] P. J. Gerrish and R. E. Lenski. The fate of competing beneficial mutations in an asexual population. 
Genetica, 102:127-144, 1998. 


37 



[22] J. B. S. Haldane. A mathematical theory of natural and artificial selection. In Mathematical Pro¬ 
ceedings of the Cambridge Philosophical Society, volume 23, page 607-615. Cambridge Univ Press, 
1927. 

[23] M. W. Hirsch. Systems of differential equations which are competitive or cooperative: lii. competing 
species. Nonlinearity, 1(1):51, 1988. 

[24] J. Hofbauer and K. Sigmund. Evolutionary game dynamics. Bulletin of the American Mathematical 
Society, 40(4):479-519, 2003. 

[25] E. C. Holmes, L. Q. Zhang, P. Simmonds, C. A. Ludlam, and A. Brown. Convergent and divergent 
sequence evolution in the surface envelope glycoprotein of human immunodeficiency virus type 1 
within a single infected patient. Proceedings of the national Academy of Sciences, 89(ll):4835-4839, 
1992. 

[26] B. Kerr, M. A. Riley, M. W. Feldman, and B. J. Bohannan. Local dispersal promotes biodiversity in 
a real-life game of rock-paper-scissors. Nature, 418(6894):171-174, 2002. 

[27] Y. Kim and H. A. Orr. Adaptation in sexuals vs. asexuals: clonal interference and the fisher-muller 
model. Genetics, 171(3):1377-1386, 2005. 

[28] B. C. Kirkup and M. A. Riley. Antibiotic-mediated antagonism leads to a bacterial game of rock- 
paper-scissors in vivo. Nature, 428(6981):412-414, 2004. 

[29] G. I. Lang, D. Botstein, and M. M. Desai. Genetic variation and the fate of beneficial mutations in 
asexual populations. Genetics, 188(3):647-661, 2011. 

[30] G. 1. Lang, D. P. Rice, M. J. Hickman, E. Sodergren, G. M. Weinstock, D. Botstein, and M. M. 
Desai. Pervasive genetic hitchhiking and clonal interference in forty evolving yeast populations. 
Nature, 500(7464):571-574, 2013. 

[31] R. Maddamsetti, R. E. Lenski, and J. E. Barrick. Adaptation, Clonal Interference, and Frequency- 
Dependent Interactions in a Long-Term Evolution Experiment with Escherichia coli. Genetics, 2015. 

[32] J. A. Metz, S. A. Geritz, G. Meszena, F. J. Jacobs, and J. Van Heerwaarden. Adaptive dynamics, a 
geometrical study of the consequences of nearly faithful reproduction. Stochastic and spatial structures 
of dynamical systems, 45:183-231, 1996. 

[33] H. J. Muller. Some genetic aspects of sex. American Naturalist, page 118-138, 1932. 

[34] H. J. Muller. The relation of recombination to mutational advance. Mutation Research/Fundamental 
and Molecular Mechanisms of Mutagenesis, l(l):2-9, 1964. 

[35] J. R. Nahum, B. N. Harding, and B. Kerr. Evolution of restraint in a structured rock-paper-scissors 
community. Proceedings of the National Academy of Sciences, 108(Supplement 2):10831-10838, 2011. 

[36] G. Neumann and S. Schuster. Gontinuous model for the rock-scissors-paper game between bacteriocin 
producing bacteria. Journal of mathematical biology, 54(6):815-846, 2007. 

[37] M. A. Nowak and K. Sigmund. Evolutionary dynamics of biological games, science, 303(5659):793- 
799, 2004. 

[38] H. A. Orr. The genetic theory of adaptation: a brief history. Nature Reviews Genetics, 6(2):119-127, 
2005. 

[39] G. E. Paquin and J. Adams. Relative fitness can decrease in evolving asexual populations of s. 
cerevisiae. 1983. 

[40] S. J. Schreiber and T. P. Killingback. Spatial heterogeneity promotes coexistence of rock-paper- 
scissors metacommunities. Theoretical population biology, 86:1-11, 2013. 


38 



[41] B. Sinervo and C. M. Lively. The rock-paper-scissors game and the evolution of alternative male 
strategies. Nature, 380(6571):240-243, 1996. 

[42] C. Smadi. An eco-evolutionary approach of adaptation and recombination in a large population of 
varying size. Stochastic Processes and their Applications, 125(5): 2054-2095, 2015. 

[43] D. R. Taylor and L. W. Aarssen. Complex competitive relationships among genotypes of three 
perennial grasses: implications for species coexistence. American Naturalist, page 305-327, 1990. 

[44] M. J. Wiser, N. Ribeck, and R. E. Lenski. Long-term dynamics of adaptation in asexual populations. 
Science, 342(6164):1364-1367, 2013. 

[45] S. Wright. Evolution in mendelian populations. Genetics, 16(2):97, 1931. 

[46] E. Zeeman and M. Zeeman. Erom local to global behavior in competitive lotka-volterra systems. 
Transactions of the American Mathematical Society, page 713-734, 2003. 

[47] M. Zeeman and P. van den Driessche. Three-dimensional competitive lotka-volterra systems with no 
periodic orbits. SIAM Journal on Applied Mathematics, 58(l):227-234, 1998. 

[48] M. L. Zeeman. Hopf bifurcations in competitive three-dimensional lotka-volterra systems. Dynamics 
and Stability of Systems, 8(3):189-216, 1993. 


39 



